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We present the Condon domain phase diagram for a silver single crystal measured in magnetic 
fields up to 28 T and temperatures down to 1.3 K. A standard ac method with a pickup coil 
system is used at low frequency for the measurements of the de Haas- van Alphen effect (dHvA). 
The transition point from the state of homogeneous magnetization to the inhomogeneous Condon 
domain state (CDS) is found as the point where a small irreversibility in the dHvA magnetization 
arises, as manifested by an extremely nonlinear response in the pickup voltage showing threshold 
character. The third harmonic content in the ac response is used to determine with high precision 
the CDS phase boundary. The experimentally determined Condon domain phase diagram is in 
good agreement with the theoretical prediction calculated by the standard Lifshitz-Kosevich (LK) 
formula. 

PACS numbers: 75.45. +j, 71.70.Di, 75.60.-d 



I. INTRODUCTION 

The formation of dia- and paramagnetic domains has 
been predicted by Condon^ to occur in non-magnetic 
pure metals by considering the collective interaction be- 
tween the electrons on Landau-quantized energy levels in 
the de Haas- van Alphen (dHvA) effect. These domains 
corresponding to an inhomogeneous magnetization are 
commonly called Condon domains. The domain forma- 
tion results from a self-consistent treatment of the oscil- 
lating dHvA magnetization M due to the orbital quan- 
tization of the electronic system in the total magnetic 
induction B = i.iq{H + M), where M — M{B) depends 
on the total induction B in an applied magnetic field 
H . Following the Pippard-Shoenberg concept of mag- 
netic interactiorf^ where the electrons experience the in- 
fluence of the magnetic field induced magnetization of all 
neighboring electrons, a thermodynamic instability arises 
when the dHvA amplitude becomes large enough, i.e., the 
differential susceptibility 
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This condition corresponds to the situation where the 
amplitude of the magnetization amplitude becomes com- 
parable to the magnetic field period of the dHvA effect. 
For sufficiently strong magnetization amplitudes, this in- 
stability condition, rewritten like fioOH/dB = 1 — x < 0, 
occurs in a certain field interval within the paramag- 
netic part of each dHvA cycle. In these field intervals, 
where fiodH/dB < 0, the induction as function of the ap- 
plied field B{H) is multi-valued, like the van der Waals 
isotherm for a real gas. The system avoids this instability 
in the same way as the real gas. For an infinite long rod- 



like sample (demagnetization factor n = 0), the induc- 
tion B undergoes a discontinuous transition between the 
two stable states with the induction Bi and B2 at a cer- 
tain applied field H, like the liquid-gas specific volumes 
change discontinuously at the equilibrium vapor pressure. 
Both stable states Bi and B2 correspond to the same 
free energy and the inductions in the instability interval 
{Bi,B2) are never realized. Figure [IJa) shows schemat- 
ically the magnetization energy Emag — 2^!^^^ ~ I^qH)^ 
for n = and the oscillating dHvA energy described by 
the LK-formula Eqsc — a cos(27rF/_B) to its simplest ap- 
proximation with a the oscillation amplitude and F the 
dHvA frequency. The sum of both energies as function of 
B is shown for three different magnetic fields in Figjljb)- 
(d). Usually, there is only one minimum in the total 
energy for a given applied magnetic field and the system 
will assume this value of B. However if the curvature 
of the Emag IS Smaller than the curvature of Eosc two 
minima coexist at an applied field Hi [Figjljc)] and the 
induction will jump discontinuously from the value Bi to 
B2 when sweeping the magnetic field through Hi. 

For a plate-like sample, oriented normal to the ap- 
plied magnetic field H [n — 1), the boundary condi- 
tion B — /io[iJ + (1 — n)M] — ^qH is required even in 
the interval Bi < fioH < B2. Therefore, the induction 
B can not change discontinuously and a homogeneous 
state is no longer possible. The plate breaks up into do- 
mains of opposite magnetization. The volume fractions 
of the domains with the respective inductions Bi and 
B2 are adjusted in a way that for the average induction 
of the sample B — ^qH is fulfillecP. The regions with 
Bi < ^qH are diamagnetic, those with B2 > iiqH are 
paramagnetic. The domain walls between the phases Bi 
and B2 run parallel to H across the plate. In contrast to 
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FIG. 1: (Color online) (a) Schematic representation of the 
magnetization energy Emag and the dHvA-energy Eosc as 
function of B for n — 0. (b)-(d) Sum of these energies for 
three applied magnetic fields Ho < Hi < H2. If the curvature 
of the parabola is smaller than the curvature of the oscillating 
energy two minima coexist at Hi. This leads to discontinu- 
ous jump of the induction when sweeping the magnetic field 
through Hi. 



magnetic domains in common ferro- and antiferromag- 
netism, Condon domains do not have their origin in the 
interaction of electrons via their spin moments, but via 
their orbital motion. 

For a sample shape with intermediate demagnetization 
factor, < n < 1, the above mentioned interval of mag- 
netic field Bi < hqH < B2 with the occurrence of the in- 
stability will be reduced compared to the plate-like sam- 
ple with the field range of this interval proportional to n. 
Therefore, samples of arbitrary shape will still show the 
nonuniform domain state with the same dia- and para- 
magnetic phases, Bi and B2, whose domain structure 
might however be more complex. 

Besides the analogy with the van der Waals gas, there 
is a close analogy of the CDS with the intermediate 
state of type-I superconductors, where the same bound- 
ary condition of a magnetic field applied to a sample of 
nonzero demagnetization factor leads to the formation of 
alternating domains in the normal and superconducting 
state. Condon domains, however, have the unique feature 
that the transition between the uniform and the inhomo- 
geneous domain state occurs periodically in subsequent 
dHvA oscillations. 

Equation [l] defines the boundary between the uniform 
and the Condon domain state. The resulting CDS phase 
diagram in the (H, T) plane can be predicted by means 



of the Lifshitz-Kosevich (LK) formula for the oscillatory 
magnetization of the dHvA signal using the Fermi sur- 
face parameters, like the curvature A" = d^A/dk"^ of the 
Fermi surface cross section A and the effective mass to* , 
and the Dingle temperature Td as a parameter for the 
impurity-scattering damping of the signaP. 

Up to now, Condon domains have been observed by 
different experimental methods: by NMR^, muon spin 
rotation (/iSR) spectroscopjEEl and, more recently, by 
local Hall probes^. All experimental observations have 
in common that two distinct inductions Bi and B2 or 
an induction splitting 5B = B2 — Bi are measured at a 
given applied field H and temperature T. However, these 
measurements yielded only a few points well inside the 
(H, T) diagram where Condon domains exist that could 
be compared with the theoretically predicted diagram. 
As a consequence, for example, the data on beryllium ob- 
tained by /xSR required new phase diagram calculations 
with a modified LK-formula for the susceptibility^. The 
exact determination of the CDS phase boundary, where 
SB approaches zero, is difficult and time-consuming with 
a difference measurement of Bi and Sji^^- 

It was shown recently that a small hysteresis occurs in 
the measured dHvA signal upon passing the CDS phase 
boundarjEl J)\xe to the irreversible magnetization, an 
extremely nonlinear response to a small modulation field 
arises in standard ac susceptibility measurements. The 
out-of-phase signal and the third harmonic of the pickup 
voltage rise steeply at the transition point to the CDS. 
The threshold character of these quantities offers there- 
fore a possibility to measure a Condon domain phase di- 
agram. One should note that the third harmonic of the 
susceptibility is commonly used as a very sensitive tool 
to detect phase boundaries also of other systems like e.g. 
the vortex-glass transition in superconductor^^. 

In this article we determine the Condon domain phase 
diagram for silver using the third harmonic of the ac sus- 
ceptibility for the detection of the nonlinear magnetic 
response. It was shown earliei^^ that detailed calcula- 
tions of the magnetoquantum oscillations in silver based 
on the LK-formula are in good agreement with experi- 
mental dHvA data up to 10 T. This is certainly due to 
the nearly spherical Fermi surface of silver. Expecting a 
good agreement with the theoretically determined CDS 
phase diagram, we applied to silver this first detailed de- 
termination of the phase diagram. 



II. EXPERIMENTAL 

The measurements were performed on a high quality 
silver single crystal of 4.1 x 2.1 x 1.0 mm^. The sample 
was cut from the same piece than the sample used for the 
direct observation of Condon domains using local Hall 
probe detectiorP. The sample preparation is described 
in detail elsewher^^^E^. The sample has a residual resis- 
tance ratio -R300 K/-R4.2 K = LBxlO"*^, measured by the 
contactless Zernov-Sharvin methotP^. The high quality 
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of the sample results in a very low Dingle temperature, 
which was estimated from standard dHvA analysis to be 
about Tj) — 0.2 K yielding an electronic mean free path 
of about 0.8 mm. 

A standard ac modulation method with a compensated 
pickup coil system was used. Both pickup coils are identi- 
cal and consist of about 400 turns. A long coil wound by 
a copper wire produced the modulation field with vari- 
able amplitude at frequencies of 20 — 200 Hz. The pickup 
voltage was simultaneously measured by two lock-in am- 
plifiers on the first and on higher harmonics. The mea- 
surements were performed in a superconducting coil up 
to 16 T as well as in a resistive coil up to 28 T at tem- 
peratures of 1.3 — 4.2 K. The long side of the sample was 
parallel to the [100]-axis of the single crystal and was 
slightly tilted (~ 5 deg) with respect to the direction of 
the applied magnetic field so that only the dHvA fre- 
quency from the "belly" orbit of 47300 T existed in the 
frequency spectrum. 

The method of nonlinear detection, we use here, is ap- 
plied to determine for the first time a CDS phase dia- 
gram over a broad range of temperatures and magnetic 
fields. Therefore, we will present carefully the technical 
details of the measurements in order to show the robust- 
ness of the phase diagram determination with respect 
to changing experimental parameters and measurement 
conditions. 



III. EVIDENCE OF HYSTERESIS IN SILVER 

The employed method to determine a Condon domain 
phase diagram is based on the appearance of hysteresis in 
the CDS which was first discovered on berylliunJiil Hys- 
teresis appears at the phase transition to the CDS and 
this results in some radical changes in the response to an 
ac modulation field. In the following, we will show that 
the characteristic nonlinear features in the ac response, 
as found in beryllium, are also observed in silver. 

In presence of hysteresis the amplitude of the suscep- 
tibility, normalized on the modulation level, depends on 
the modulation amplitude. This is expected to occur 
when the modulation level is of the order of the hysteresis 
loop width. The schematic representation of a hysteresis 
loop in Fig. [2] explains this nonlinear response to an ac 
field modulation. As a result, after the transition to the 
CDS, the positive (paramagnetic) part of the susceptibil- 
ity turns out to be reduced. From a comparison of two 
normalized susceptibilities, one measured with high and 
the other with low modulation level, we can in principle 
find where the amplitude reduction starts and thereby 
the transition point to the CDS. 

Figure [2] shows as well that the response to a sinusoidal 
field modulation becomes window shaped and is slightly 
shifted in phase with respect to the input. Therefore, 
both the third harmonic and the out-of-phase signal of 
the pickup voltage increase steeply when the CDS phase 
boundary is crossecff^'. This threshold behavior offers 
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FIG. 2: (Color online) Schematic representation of two minor 
hysteresis loops. The ac response upon field modulation is 
window-Hke and slightly shifted in phase. The normalized 
response (apparent susceptibility) decreases sharply when the 
modulation amplitude becomes of the order of the hysteresis 
loop width. 



a simple way to determine the transition point of the 
CDS. The major advantage of third harmonic and out-of- 
phase part measurements is that only one magnetic field 
or temperature sweep through the transition is needed. 

Figures [3] and |4] show the above discussed nonlinear 
features in the pickup signal at constant temperature 
T = 2.7 K measured in the superconducting coil. Fig- 
ure [3]j a) shows two traces of the normalized pickup volt- 
age, i.e. the susceptibility, obtained in the same con- 
ditions with 1.0 G and 0.2 G modulation amplitude at 
160 Hz modulation frequency. In principle, both mod- 
ulation levels are small enough compared to the dHvA 
period of about 20 G at 10 T that identical traces are ex- 
pected for the susceptibility. The expanded view around 
8 T, which is outside the CDS, shows that the normal- 
ized signals are indeed identical. For higher fields, on the 
other hand, the upper part of the susceptibility wave- 
form, measured with the smaller modulation level, is re- 
duced. The expanded view around 13.5 T shows that the 
signals are identical except for the positive part of the 
oscillation. This implies that at this part of the dHvA 
oscillation the magnetization is irreversible and there is 
a small hysteresis loop. The width of the hysteresis loop 
is of the order of 0.2 G. 

A similar decrease of the normalized response was ob- 
served earlieJi^ on silver at low temperatures. Unfortu- 
nately, because of the absence of the experimental pa- 
rameters, this study can be compared only qualitatively 
with our data. 

The magnetic field where the normalized pickup volt- 
ages start to differ between low and high modulation level 
is marked approximately by an arrow in Fig. |3][a). We 
obtain for the critical magnetic field fioHci = 10.0 T. 

Figure [sj^b) shows the behavior of the third harmonic 
which was simultaneously measured with the first har- 
monic response in Fig. p[a) for 0.2 G modulation am- 
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FIG. 3: (Color online) (a) Pickup voltage normalized on the 
modulation level for low and high modulation level. Up to 
about 10 T the response is linear with respect to the modula- 
tion level, (b) Third harmonic of the pickup voltage measured 
at 0.2 G modulation amplitude showing that starting from 
9.5 T the harmonic content in the response increases steeply. 
Lower part of the figure shows respective zooms. Both data 
measured at 2.7 K. 



plitudc. For magnetic fields lower than the critical field 
there is only noise. At the transition to the CDS hystere- 
sis arises and the third harmonic increases very steeply. 
This is nicely seen in the respective expanded views. 
The critical field of the CDS phase boundary can be 
obtained as the intersection of the two straight lines 
shown in Fig. |3jb). Here, the critical field is found as 
tJ-oHc2 = 9.6 T. 

The amplitude of the third harmonic is expected to go 
to zero in each diamagnetic part of the dHvA period be- 
cause the sample magnetization is here homogeneous and 
without hysteresis. The presented behavior in Fig. [sf^b) 
does not go to zero exactly which is certainly the result 
of a small rectification effect or, what is the same, the 
result of phase smearing of the oscillation signal. The 
homogeneity of the coil is about 30 ppm in a sphere with 
1 cm diameter which may result in a field inhomogeneity 
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FIG. 4: (Color online) (a) Phase angle of the third harmonic 
showing clearly the transition between noise outside the CDS 
and a fixed phase in the CDS. (b) The out-of-phase part of the 
first harmonic response changes due to the arising hysteresis. 
Both data measured at 2.7 K and 0.2 G modulation level. 



of about 1 G in the sample volume at 10 T. Therefore, 
the transition to the CDS does not occur simultaneously 
in the whole sample. This effect will be much bigger in 
a resistive coil where the homogeneity is 20 times worse. 
However, we will see below, that the third harmonic rec- 
tification does not affect the determination of the critical 
field of the CDS. 

In Fig. Qa) the phase angle of the third harmonic is 
shown for the same conditions like in Fig. [3] For mag- 
netic fields where the amplitude of the third harmonic is 
below the noise level its phase angle is not determined. 
Therefore, the phase varies between —180 to -1-180 de- 
gree. With the appearance of a third harmonic signal at 
the transition to the CDS the phase becomes finite. This 
passage has a threshold character as well. The arrow in 
Fig. |4|^a) shows the position of the threshold which yields 
the critical field fioHc3 — 9.7 T. 

The behavior of the out-of-phase part of the first har- 
monic response, shown in Fig. |4|^b), offers another pos- 
sibility to determine the critical field. In the uniform 
state, without domains, the imaginary part is small and 
varies smoothly especially at low magnetic field due to 
the magnetoresistance and changing eddy currents. After 
the transition to the CDS the out-of-phase signal changes 
rapidly. The transition point can be found as the inter- 
section of two lines, as it is shown in Fig. |4|b). Here, we 
obtain for the critical field HqHca = 9.4 T. 

A comparison of the values -ffci...c4 for the transition 
field at 2.7 K shows that they are very close. We note that 



the critical field of about 10 T agrees roughly with the 
phase boundary found in the Hall probe experiment^. 
All above presented methods could, in principle, be used 
to determine the phase boundary of the CDS. The first 
method (fig.jsjja)) requires at least two field sweeps. Mea- 
surements of the out-of-phase part (fig. |4|^b)) are not 
precise, due to the high conductivity of silver and the 
resulting eddy currents (The situation might be different 
in a less conducting metal). Therefore, for silver the third 
harmonic measurements to determine the CDS phase dia- 
gram are preferred. Moreover, we will see below that the 
obtained values with the third harmonic for the phase 
boundary {H, T) do not depend drastically on the fre- 
quency and amplitude of the field modulation which of- 
fers the possibility to measure in noisier conditions. The 
found scattering in the values of the transition fields ob- 
tained from the different methods gives an uncertainty of 
about ±0.5 T in the transition fields. 



IV. PHASE DIAGRAM 
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FIG. 5: (Color online) Third harmonic response for two mod- 
ulation levels measured in a superconducting coil at a modu- 
lation frequency of 160 Hz at 2.7 K. The same critical field is 
found from the steep increase. 



Because of the increased noise level in the water-cooled 
resistive magnets, we needed to increase the signal-to- 
noise ratio by using rather high modulation frequencies 
wl60 Hz and higher modulation amplitudes, 1.0 G and 
more. In the following we check whether the modulation 
frequency and amplitude can be varied without chang- 
ing the value of the critical field deduced from the third 
harmonic response. 

Modulation amplitudes of the order of the width of the 
hysteresis loop are required to resolve the amplitude re- 
duction in the normalized pickup voltage in Fig.jsjja). For 
the third harmonic signal, as shown in Fig. [2] the non- 
linear features persist up to high modulation amplitudes. 
Figure [5] shows traces with 0.2 and 1.0 G modulation am- 
plitude. For 1.0 G modulation amplitude there is a very 
small third harmonic signal before the transition to the 
CDS takes place. This small contribution to the third 
harmonic is due to the nonlinearity of the dHvA effect 
itself'^. Nevertheless, the position of the sharp increase 
remains unchanged. 

Figure [6] shows that increasing the modulation ampli- 
tude up to 10 G and varying the modulation frequency by 
a factor four between 40 Hz and 160 Hz does not change 
the position of the critical field, either. The results pre- 
sented here were obtained in the resistive magnet. The 
measurements were made at low temperatures in order to 
compare them with data obtained in the superconducting 
magnet. 

All results for the CDS transition points obtained in 
the superconducting and the resistive magnets are pre- 
sented in Fig. [7] The critical fields for each temperature 
are found as the field where the third harmonic response 
starts to arise like in Fig. [6] One should note that near 
the flat maximum of the phase diagram T-sweep mea- 
surements would be in principle better. The solid line in 
Fig.[7]is the CDS boundary calculated for silver using the 
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FIG. 6: Third harmonic response measured after averaging 
over the oscillations in a resistive coil for different modulation 
amplitudes (a) and modulation frequencies (b).The critical 
field values obtained at a given temperature are independent 
on modulation frequency and amplitude. 



LK-formula for the susceptibility criterion x = 1 with a 
Dingle temperature of Td = 0.2 K for our sampld^^l^. 
For comparison, theoretical diagrams for Td = 0.1 K and 
0.8 K are shown in Fig. [Tj 

A good agreement of the phase diagram predictions 
based on the LK-formula with our data can be seen for 
a Dingle temperature Td = 0.2 K. Data points obtained 
in superconducting and resistive magnets overlap which 
supports that the different measurement conditions did 
not affect the precise determination of the phase bound- 
ary 
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FIG. 7: (Color online) Phase diagram in the {H, T) plane 
for silver. Experimental points from the superconducting and 
resistive magnet. The solid line is the CDS boundary calcu- 
lated by the LK-formula for Td = 0.2 Kl^^. The dashed and 
the dotted lines correspond to Td = 0.8 K and Td = 0.1 K, 
respectively. 



domain state. This substantiates that hysteresis is likely 
to occur in all pure metals that exhibit Condon domains. 
The hysteresis leads to with threshold character arising 
extremely nonlinear response to a modulation field. In 
particular, the third harmonic response to a modulation 
field increases sharply upon entering into the Condon do- 
main state. This offered the possibility to determine eas- 
ily the CDS phase boundary with high accuracy. 

The critical fields obtained from the third harmonic of 
the pickup signal of the ac modulation technique, turned 
out to be independent on changes of the modulation fre- 
quency and amplitude. Due to this independence this 
method could be used with higher modulation frequen- 
cies in pulse magnetic fields. 

Very good agreement of the CDS phase diagram is 
found with calculations of the dHvA signal based on the 
LK-theory. This agreement shows that the LK-formula 
describes well the field dependent magnetization of the 
nearly spherical Fermi surface of silver. Furthermore, 
the agreement demonstrates that the described method 
is correct for the determination of the CDS phase dia- 
gram. 



CONCLUSION 
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